Why three-body physics does not solve the proton radius puzzle 
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The possible involvement of weakly bound three-body systems in the muonic hydrogen spec- 
troscopy experiment [ij , which could resolve the current discrepancy between determinations of the 
proton radius, is investigated. Using variational calculations with complex coordinate rotation, it 
is shown that the p^e ion, which was recently proposed as a possible candidate 0], has no resonant 
states in the energy region of interest. QED level shifts are included phenomenologically by includ- 
ing a Yukawa potential in the three-body Coulomb Hamiltonian before diagonalization. It is also 
shown that the ppfj, molecular ion cannot play any role in the observed line. 

PACS numbers: 36.10.Ee 31.15.ac 31.15.xt 



Introduction. The recent measurement of the Lamb 
shift in muonic hydrogen [l| , which resulted in a new de- 
termination of the proton rms charge radius, has given 
rise to an abundant literature attempting to resolve the 4- 
percent discrepancy with its current CODATA value Q, 
mainly obtained from electron-proton (e-p) scattering 
and hydrogen (H) spectroscopv. The discrepancy with 
e-p scattering measurements may be less severe than 
the initially announced 5 standard deviations due to a 
possible underestimation of error bars in e-p data anal- 
ysis [sl. However, no agreement with H spectroscopic 
data j^] could be obtained unless some unknown system- 
atic shift was discovered in all these experiments. An 
attempt to solve the discrepancy by re-evaluating the 
third Zemach moment of the proton Q was shown to 
go against e-p scattering data 8] . Reconsideration of the 
assumptions on form factors used to extract the proton 
radius only led to a slight increase of the uncertainty . 
Much effort has been focused on checking and improv- 
ing the Lamb shift calculations used to extract the pro- 
ton radius (see e.g. 0, [13] and references therein, as well 
as [HI). In particular, the polarizability of the proton 
is a subject of hot debates [l^], but no consensus has 
been reached yet regarding claims that this contribution 
may be able to solve the proton size puzzle. The muonic 
hydrogen result also triggered a search for new physics, 
which is however highly constrained by many low-energy 
data 

So far, no possible experimental artifacts have been 
proposed, with the notable exception of [I]- In that ref- 
erence, Jentschura suggested that ^j,p atoms formed in 
the 2S state might bind an electron present in the gas 
target to form a (pfie)^ ion. Since the muon orbit is 
close to the proton, the proton charge is shielded from 
the outer electron, but attractive dipole interactions with 
the /ip(2S) core could result in a quasibound (or reso- 
nant) state slightly below the 2S threshold. Jentschura 
gave an order-of-magnitude estimate of the associated 



frequency shift, which was indeed compatible with the 



observed discrepancy AE^. 
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tween experiment and theory. Among the many proposed 
solutions to the proton radius puzzle, this appears to be 
one of the last unexplored paths, and definitely deserves 
further study. The purpose of Section U is to give a con- 
clusive -eventually negative- answer on the viability of 
the Pile hypothesis. 

Other candidates to form such quasibound states with 
a /ip atom can also be considered. While the possibility of 
forming (p/i/i)^ ions can be safely ignored -if only because 
of the low number of niuons in the experiment chamber, 
it is known that (pp/i)"*" molecules can be formed in col- 
lisions of excited /ip(2S) atoms with hydrogen molecules 
(see e.g. [ij] and references therein). In this process, 
the favored molecular states lie ~ 10 — 20 eV below the 
n = 2 dissociation threshold, and their Coulomb lifetime 
is much too short to account for the observed transi- 
tions [l^ [l^. However, it is worth reconsidering the 
possible role of more weakly bound states which can have 
long lifetimes. In Section [Hi we discuss the relevant por- 
tion of the spectrum of pp^i and show that none of the 
existing states has the required characteristics to match 
the experimental data. 



I. THE p/xe HYPOTHESIS 

A. Coulomb Schrodinger Hamiltonian 

A first insight into the physics of the p^ie system is pro- 
vided by the adiabatic approximation: the "slow" outer 
electron experiences the attractive dipole potential cre- 
ated by the p/i(n = 2) core. Neglecting (in a first step) 
the 2S-2P Lamb shift, the n = 2 degenerate level has 
a nonzero static dipole A = Sme/rrir in atomic units, 
where rrir is the pfi reduced mass. The long-range behav- 
ior of the potential curve is then —A/R^, and it has been 
known for a long time [l7j that in this kind of potential, 
a particle of mass m has an infinite series of bound states 
with exponentially decreasing binding energies, provided 
A > Ac = me/^m. This translates into the condition 
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m > mr/24 ~ 7.74 irie- Since the p^e system is below 
that value, only a finite number of bound states can ex- 
ist. This is a first indication that due to the weakness of 
the attractive dipole potential, the existence of resonant 
states is not assured. 

For a more detailed investigation, we consider the full- 
three body dynamics and use a numerical diagonalization 
of the complex rotated 18] Schrodinger Hamiltonian in a 
variational basis set to obtain the energies and widths of 
resonances. Having in mind that our goal is to determine 
whether resonant states do exist, we start with the p^^ 
negative ion, where resonances below the n = 2 thresh- 
old are known to exist, and then follow the position of 
the lowest resonant state (with ro-vibrational quantum 
numbers v — L — 0) while decreasing the mass m of 
the third particle -which we denote x in the following- in 
small steps, from the muon down to the electron mass. 
We used the perimetric coordinates (denoted x, y, z) and 
Sturmian wave functions 



xL„^ [ax) L„„ (/3y) L„^ (/3z) e-("-+/'?^+/'-)/2 ^ ' 

where nx,ny, Uz are non- negative integers and L„ the La- 
guerre polynomials. The main idea of this method (Tgj 
is to take advantage of the dynamical symmetries of the 
three body Coulomb problem; in this basis set the Hamil- 
tonian has strict coupling rules, leading to a sparse-band 
matrix which can be diagonalized in fewer operations 
with respect to a "full" matrix, maintaining better nu- 
merical stability 0, . This choice is well suited for 
the present study, in particular since the small number 
of variational parameters (the two length scales and 
makes it easy to readjust them for each value of to. 

Our results are summarized in Fig. [1] Numerical ac- 
curacy is better than 10~^ eV, so that even the small- 
est width, lying in the 10~* eV range, is obtained with 
two significant digits. The point at the extreme right 
corresponds to the pfifj, negative ion; its binding energy 
Et = 120.4716427 eV and width F = 7.146063657 eV 
(obtained with the 2010 CODATA values of the proton 
and muon masses ,3]), improve on the values reported 
in^, Eb = 120.39 eV and F = 7.12 eV. As expected, 
when the mass of the third particle decreases, its orbit 
around the pfx{n = 2) core becomes larger and its binding 
energy decreases. 

Coulomb widths are shown mainly for illustration since 
the focus is on the existence of quasibound states. They 
decrease together with Ef, with decreasing to, because 
the system comes closer and closer to the adiabatic limit 
where the resonant state tends to a true bound state. 
The sharp drop between to = 48 to^ and to = 47 me 
is due to a crossing of the pii{n = 2) and px{n — 1) 
thresholds which occurs at to = TOj,to^/(4top -|- 3TOp) ~ 
47.67 TOe. Below that value, the px{n — 1) threshold is 
above p/i(n = 2), so that one of the channels for Coulomb 
deexcitation becomes forbidden, and the width decreases. 

No calculations were performed below to = 32TOe, 
where larger computer resources would be needed to con- 
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FIG. 1: Binding energy (top) Et — E{pfi)„=2 — E and 
Coulomb width (bottom) of the lowest ^S'^ resonance of the 
p'^H'x' system, as a function of the mass m of the third 
particle x. Squares: pure Coulomb Hamiltonian, circles: 
Coulomb Hamiltonian with an added Yukawa potential ad- 
justed to reproduce the 2S-2P Lamb shift. For the latter, 
the binding energy is with respect to the shifted 2S thresh- 
old. On the top graph, the horizontal line recalls the order 
of magnitude of the discrepancy between theory and experi- 
ment in the proton radius experiment, while the vertical line 
at m/nie — 7.74 signals the minimal value of m to retain an in- 
finite series of resonant states in the adiabatic approximation. 
On the bottom graph, the horizontal line shows the value of 
the width corresponding to a Coulomb lifetime Tc — 1 fj,s, 
and the vertical one at m/nie = 47.67 shows the crossing 
point between the pfJ,{n — 2) and px{n = 1) thresholds. 



verge the resonance width with sufficient accuracy. How- 
ever, the results of Fig. [l]are amply sufficient to demon- 
strate that for TO = TOe the binding energy of a resonant 
state would be too small by many orders of magnitude to 
resolve the 0.31 meV discrepancy in the muonic hydrogen 
experiment. 



B. Inclusion of QED effects 

An important ingredient is still missing in our treat- 
ment, that is the Lamb shift which lifts the degeneracy 
of the n = 2 level, leading to a zero static dipole. The 
dipole potential in —A/R'^ is thus replaced by an induced 
dipole potential decaying like —ao/2R'^ where ao is the 
polarizability of the 2S (or 2P) state. As a consequence. 
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when resonances come closer and closer to the threshold, 
they experience a more and more weakened attractive po- 
tential. The infinite dipole series obtained for po- 
tentials becomes truncated since a potential only 
admits a finite number of bound states [l^, [2^. It is 
therefore clear that the QED shift of energy levels does 
not put into question the above conclusion, since it weak- 
ens the attractive potential and can only favor the dis- 
appearance of weakly bormd resonant states. 

In order to show this more quantitatively, we repre- 
sent QED corrections by an effective potential added to 
the three-body Coulomb Hamiltonian before diagonal- 
ization (the perturbative approach is not valid here due 
to the mixing of 2s and 2p orbitals). The most nat- 
ural choice would be the Uehling vacuum polarization 
potential which accounts for more than 99% of the /ip 
Lamb shift. This approach was used to obtain the en- 
ergy shifts of dt^, resonances using the stabilization tech- 
nique with a real scaling parameter [22]. However, that 
method would not yield precise enough results here be- 
cause the states under consideration have a large width 
that is not negligible with respect to their binding en- 
ergy. In principle, the complex rotation method allows 
for higher precision, but numerical experimentation sug- 
gests that the Uehling potential is not dilation analytic 
(see e.g. [23| for a definition of this term) and thus cannot 
be easily combined with the complex rotation method. 
For this reason, we "mimic" the proton- muon Uehling 
potential by a Yukawa potential Veff{r) = — Ae~"''/r. 
First, we determined the exact shifts of the 2S and 2P 
levels in jip induced by the Uehling potential by diag- 
onalization of the two-body Coulomb-|-Uehling Hamil- 
tonian: AE{2P) = -14.5789 meV and AE{2S) = 
-219.7374 meV The deduced contribution to the 

Lamb shift, AE{2P) - AE{2S) = 205.1584 meV is in 
exact agreement with the result of 2]. Then, the pa- 
rameters a = 492.86 and A = 2.050982 10"^ (in atomic 
units) were adjusted so that the diagonalization of the 
two-body Coulomb+Yukawa Hamiltonian yield the same 
exact shifts. 

Going back to the three-body problem, the perimet- 
ric coordinates method loses its point here, because the 
added potential, having no selection rules, is represented 
by a "full" matrix. We thus used another variational 
expansion with Hylleraas coordinates ri,r2,ri2 and ex- 
ponential basis functions 

xT\n,r2,n2) = Re(Im)[e-""'---^"'-^-^"'-^] (2) 

with pseudorandom complex exponents a„,/3„,7„ ^2^. 
In order to test the validity of the Yukawa approxima- 
tion, we computed the energy of the lowest resonances of 
dtfi below the n = 2 threshold. Agreement with calcu- 
lations using the exact Uehling potential [22| is obtained 
at the meV level (see Table IJ) . The approximation is ex- 
pected to work even better for more weakly bound res- 
onances, since the outer particle's orbit lies farther from 
the np{n — 2) core and will be less sensitive to the short- 
range details of the proton-muon interaction potential. 
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TABLE I: Binding energies of resonant states of dt/i (top) and 
p^x (bottom) below the /ip(2S) threshold. Eb is obtained with 
the pure Coulomb Hamiltonian, i?^ with an added Yukawa 
potential (see text), and E^ (for dt^ only) with the exact 
Uehling potential. Resonance positions are given with respect 
to the shifted 2S atomic threshold. For p/ix, the corresponding 
widths F, F^ are also given. 

The main approximation in our approach is the neglect 
of fine and hyperfine structure (e.g. the 2Si/2 hyperfine 
splitting is 22.8 meV). This is not expected to strongly 
modify our results, since the polarizabilities of individual 
hyperfine levels differ from the structureless case by less 
than 2 percent. We estimate that QED-induced changes 
in binding energies and widths are obtained with a rela- 
tive accuracy of a few percent. 

Results are shown in Table|T]and plotted in Fig. [T] (open 
circles). The drop of binding energy ranges from about 
62 meV for m = 80 nif, down to a few meV for to = 50 TOg. 
This is consistent with the picture of a resonant state as 
a jip atom in a 2S orbital interacting with the field of the 
distant x particle. Indeed, the 2S-2P mixing decreases 
as X becomes more loosely bound so that the shift of 
the resonance tends to that of the 2S level (220 meV). 
Already for m — A&nrie, the binding energy is in the 0.1 
meV range. These results confirm that the p/ie system 
can have no resonance that would be sufficiently shifted 
with respect to the 2S atomic level to explain even a small 
fraction of the 0.3 meV discrepancy. 



II. MOLECULAR STATES 

As mentioned in the Introduction, metastable states 
of the pp^ molecule are known to play a role in the 
muonic hydrogen deexcitation cascade [3, 1131 • How- 
ever, their possible involvement in the observed transi- 
tion signal is highly constrained by experimental data [3| . 
Firstly, the observed 18 GHz-wide transition cannot be 
linked to a photo-association (or photo-dissociation) pro- 
cess, which would lead to a much broader spectrum. Sec- 
ondly, two hyperfine components were observed, and the 
second one [28| was found at the predicted frequency 
with the proton radius value deduced from the first one 
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i.e. the measured line has a hyperfine structure that 
closely matches that of the 2Si/2 2P3/2 atomic transi- 
tion. Thus the only possible scenario involving molecular 
states is that of a transition between resonant states, a 
first one lying close to the /xp(2Si/2) threshold and a sec- 
ond one close to the /ip(2P3/2) threshold. A prerequisite 
to a resolution of the proton radius puzzle would be to 
find such a pair of states, connected by a dipole-allowed 
transition, with such binding energies that the transition 
would be blue-shifted by a fraction of meV with respect 
to the 2S1/2-2P3/2 atomic transition. One strong point 
of this hypothesis is that the lower state would have long 
enough Coulomb and radiative lifetimes (since its wave 
function mainly involves the 2S atomic orbital with ex- 
tremely small admixture with 2P) to survive the ^ 1/is 
delay before the spectroscopy laser pulse; similarly, the 
upper state wave function being mainly of 2P character, 
the transition width could be close to the observed one. 




Internuclear distance (a ) Internuclear distance (a ) 



FIG. 2: (a) Binding BO electronic energy curves below the 
n = 2 threshold. Distances are in units of the muonic Bohr 
radius ~ ao/207, and energies are with respect to the 
threshold, (b) Enlarged view of the region of large internu- 
clear distances. The thin black lines are the continuation of 
the curves in (a), and converge to the unshifted n = 2 thresh- 
old. The vertical scale on the left (right) correspond to the a 
(tt) curves respectively. The a curves are indistinguishable at 
such large distances; their convergence is slower because they 
behave like instead of 1/-R*. The thick red lines are 

obtained by adding a Yukawa potential Y representing the 
proton-muon Uehling potential. The 2p7r„ curve converges to 
the shifted 2P threshold shown on the right vertical axis. The 
3dag and 4/cti, curves converge to the shifted 2S threshold 
shown on the left axis. Since BO curves are drawn for infinite 
nuclear masses, for consistency the vacuum polarization shifts 
and associated parameters of Y were calculated for a fip atom 
with infinite proton mass: AE{2S) = —281.5896 meV and 
AE{2P) = -27.5666 meV yielding an exponent a = 443.7624 
and coefficient A — 1.919722 10~^ for the Yukawa potential 
(see Sect. IB}. 

There are three binding Born-Oppenheimer (BO) elec- 
tronic curves below the n — 2 threshold (see Fig [5^) . 
The 3dag and 4/f7„ curves have a, 1/B? long-range be- 
havior (dipole potential) and support infinite series of 



resonant states of symmetries ^P°, ^D'^... and ^S^ , 
^D"... respectively [ll,[il,[2|. The 2p7r„ curve de- 
cays like (induced dipole potential) and supports 
a finite number of states with symmetries , ^P'^ ^ ^D°, 
^D"... (30l43^. The Lamb shift splits the 71 = 2 threshold, 
and the key point is that both a curves converge to the 
fip{2S) threshold, while only the 2pT:u curve converges 
to the 2P threshold. This is clearly apparent in Fig[5}3, 
which shows how the BO energies are modified by the 
Lamb shift at large internuclear distances. We obtained 
these curves by calculating the energy of an electron in 
the field of the nuclei for a grid of values of R, with a 
Yukawa potential added to the Coulomb potential in or- 
der to represent the Uehling potential, as described in 
Sect. ITB] The electronic wave function was expanded on 
a variational set of exponential basis functions with pseu- 
dorandom exponents [33| . 

These considerations show that the upper level of the 
proposed molecular transition could only be one of the 
levels supported by the 2p7r„ curve. These states have 
already been investigated in detail [30|-|32|: one res- 
onant state (with a binding energy E}, = 13.57 eV) and 
one "anomalous parity" ^P*^ bound state (i?{, = 13.54 
eV ) ha ve been evidenced, and it was shown convincingly 
in that no other state can exist fs^]. We can thus 
conclude that the hypothesis of a molecular transition 
fails due to the lack of a suitable upper state. 

Conclusion. We have refuted all reasonable hypothe- 
ses aiming to resolve the " proton radius puzzle" with the 
help of three-body physics. On the one hand, it was 
shown that the p/ie atom has no resonant states below 
the iJ,p{n = 2) threshold; on the other hand, that the 
spectrum of the ppjj, molecular ion leaves no room for 
a transition lying close to the 2S-2P atomic line, since 
there are no long-lived resonant states in the vicinity 
of the 2P threshold. However, the search for possible 
experimental artifacts in the muonic hydrogen experi- 
ment might not be over yet; as noted in formation 
and deexcitation of muonic hydrogen are very compli- 
cated processes [sS] in which other many-body bound or 
quasibound states could play a role. Another result of 
this work is the demonstration of a convenient method 
to study the effect of QED level shifts on the position 
and width of few-body resonances in the framework of 
complex-rotated variational methods, using an adjusted 
Yukawa potential. 
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